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Abstract 

Employing the heavy quark pole mass and static potential that were obtained from resummation 
of the renormalon-caused large order behavior to all orders, and matching the resummed static 
potential to the long distance Cornell potential we obtain a new estimate of the bottom quark mass 
'nI" ■ from T(1S) system. We point out that the pure perturbative QCD approach to T(1S) spectrum 

^\ . should have a sizable systematic error of about 100 MeV, due to the absence of the long distance 

CN ', confining potential, and that the incorporation of the latter is essential for accurate determination 

Ch I of the ground state energy. We also recalculate the nonperturbative Leutwyler-Voloshin effect 

using a resummed octet potential matched to a confining linear potential, and obtain a result that 
is substantially smaller than the estimates based on Coulombic potentials. The hyperfine splitting 
CN '. is estimated to be 50 ± 8 MeV at a^iM,,) = 0.118. 
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I. INTRODUCTION 

The heavy quarkonium is a quantum chromodynamic (QCD) system that essentially 
depends on the two fundamental parameters, the heavy quark mass and strong coupling 
constant. The precise experimental determination of the spectra of heavy quarkonium such 
as bottomonium provides a good opportunity to test QCD as well as determine the funda- 
mental parameters. 

Traditionally, the heavy quarkonium was investigated primarily within the phenomeno- 
logical potential models that proved very successful. Nevertheless, the ambiguities in the 
potential models forbade a precise connection to QCD, and thus an accurate determination 
of the QCD parameters in this approach was not feasible. The recent approaches to quarko- 
nium system is based on effective theories such as the potential NRQCD (pNRQCD) jy. In 
pNRQCD Lagrangian the heavy quark mass, which appears in its most natural form as a 
pole mass, and the interquark potentials are matching coefficients, which can be calculated, 
formally, in perturbative QCD (pQCD). 

For the bottomonium the first energy eigenstates are expected to be largely determined 
by the interquark potentials at short distances which can be calculated in pQCD. If the 
nonperturbative effects are indeed small then these first states are essentially perturbative, 
and can open a window for testing pQCD as well as determining the QCD parameters. 
Crucial for the success of this approach is that the nonperturbative effects be small and the 
potentials at short distances as well as the pole mass be calculable reliably within pQCD. 
It is, however, well known that because of the infrared (IR) renormalons the perturbative 
expansions for the pole mass as well as the static potential, which are known to next- 
next leading order (NNLO), show bad convergence even at scales where pQCD should be 
applicable. 

Besides the apparent convergence problem the IR renormalons pose a more fundamental 
problem of defining the infrared-sensitive quantities. Let us focus first on the pole mass for 
the moment, but the discussion in the following applies equally well to the static potential. 
The renormalon divergence in the perturbative expansion of the pole mass causes an intrinsic 
uncertainty of 0{Kqcd)-, and this results in unavoidable ambiguity in defining the pole mass. 
Nevertheless, a pole mass can be well defined by Borel resummation. In Borel resummation 
the renormalon ambiguity appears as an ambiguity in taking the integration contour over 
the renormalon singularity on the Borel plane. Despite this ambiguity one can define a pole 
mass simply by taking the principal value prescription for the Borel integral. We call the 
pole mass defined this way the BR mass J2, |3[. The BR mass, by definition, has exactly 
the same perturbative expansion as the pole mass in pQCD, and thus the coefficients suffer 
from the same renormalon divergence, but nonetheless it is well-defined. Obviously, the 
BR mass is not a short- distance mass, and in this respect it is unique among the various 
heavy quark mass definitions. The ambiguity in defining the pole mass cancel out in physical 
obeservables, for instance in heavy quarkonium Hamiltonian the ambiguities in the pole mass 
and the static potential cancel each other, leaving the Hamiltonian free from renormalon 
divergence y, H |6|. It is noted that we focus, throughout the paper, exclusively on the 
first IR renormalon since it is overwhelmly dominant in the divergence of the perturbative 
coefficients of the pole mass as well as the static potential. 

Although the Borel resummation formally resolves the renormalon problem, practicing 
it is another matter since the calculation of the Borel integral to an accuracy better than 
the renormalon ambiguity demands the Borel transform be known precisely in the domain 



containing the origin as well as the first IR renormalon singularity. Although the behav- 
ior of the Borel transform about the origin can be obtained by the ordinary perturbative 
expansions, obtaining its behavior from the origin to the renormalon singularity with the 
known first terms of perturbation is nontrivial, since it requires accurate knowledge on the 
renormalon-caused large order behavior. This difficulty of Borel resummation led in the 
past to the abandonment of an independent notion of infrared sensitive quantities such as 
the pole mass, and to the adoption of the short- distance quantities instead [3]. 

Since the difficulty of achieving a manageable definition of an infrared sensitive quantity 
arises from the difficulty of describing the Borel transform about the renormalon singularity, 
it is worthwhile to look at this problem more carefully. As is well known, the exact nature of 
the renormalon singularity can be determined by renormalization group (RG) equation, and 
the only missing information about the renormalon singularity is the residue that determines 
the overall normalization constant of the large order behavior. Although the residue cannot 
be calculated exactly, there is a perturbation scheme in which the residue can be expressed 
as a convergent series that can be calculated using the ordinary perturbative expansions 
[8|,|9|. The scheme, when applied to QCD observables, yields residues that show surprisingly 
good convergence even with the first terms of the perturbative expansions, which indicates 
the residues can be calculated fairly well from the low order perturbative calculations. It 
was observed recently that this method works particularly well for the pole mass and the 
static potential because of the overwhelming dominance of the first IR renormalon in these 
quantities 10|, |ll|, |l2 . 



With the residue known approximately, one can obtain a good description of the Borel 
transform in the region that contains the origin and the renormalon singularity through an 
interpolation of the perturbative expansions about the origin and the functional behavior 
about the renornalom singularity. A particular prescription of the interpolation called bilocal 
expansion was proposed in 13|, which incorporates the correct nature of the renormalon 



singularity and gives, by definition, exactly the same perturbative expansion about the 
origin as the true Borel transform. The Borel resummation employing this interpolating 
Borel transform then sums the renormalon-caused large order behavior to all orders, giving 
an all-order definition for the infrared sensitive quantity. The pole mass and static potential 
computed in this scheme show remarkable convergence under the bilocal expansion, and in 
the case of the latter an excellent agreement was obtained with the lattice potential [l3|. 
With the Borel summation technique we thus have a well-defined pole mass that has exactly 
the same perturbative expansion as the pole mass in pQCD, but nonetheless converges 
rapidly under the bilocal expansion. 

In this paper we apply the Borel resummed pole mass and the static potential to T(1S) 
system to determine the bottom quark mass and the hyperfine splitting. Our approach is 
unique in that it is not based on short-distance quantities. An useful feature of employing 
the all-order quantities is that a natural scale separation occurs in a multi-scale system such 
as the quarkonium. In the schemes that employ short- distance quantities the large and small 
scales mix invariably in the perturbative calculations, causing large perturbative coefficients 
that render the calculations less reliable. Since in our scheme the pole mass and the potential 
can be resummed at scales that are optimal for each quantities, for example, the heavy 
quark mass and the inverse of the interquark distance, respectively, the Hamiltonian in our 
approach does not suffer from the scale mixing problem. This is a feature that differentiates 
our scheme from those based on short-distance quantities. There are other advantages, for 
instance, having a well-defined, stand-alone potential makes it easy to merge the resummed 



potential at short distance to the phenomenological long distance confining potential. 

The paper is organized as follows. In Sec. II, we define the Hamiltonian in BR scheme that 
is obtained, at short and intermediate distances, rigorously from pQCD plus resummation 
and at intermediate and long distance by matching the long distance Cornell potential 
to the BR potential. In Sec. Ill, we show that the incorporation of the long distance 
confining potential is important for T(1S) and the computations based on pure pQCD 
should contain sizable systematic error. In Sec. IV, we recalculate the Leutwyler-Voloshin 
type nonperturbative effects employing the resummed singlet and octet potentials and show 
that they are sufficiently small enough to be negligible for the IS state. In Sec. V we give 
our estimate of the MS bottom quark mass and hyperfine splitting, and in Sec. VI give a 
summary. 

II. HAMILTONIAN 

The pNRQCD Lagrangian for heavy quarkonium contains singlet as well as octet fields 
and ultrasoft gluons which represent the dynamical degrees of freedom at scales below the 
typical heavy-quark momentum. Ignoring for the moment the interactions between the 
singlet and octet fields, the pNRQCD Hamiltonian for singlet quarkonium in on-shell scheme 
is given by IJ 



H = Ho + H\ (1) 

where 

H, = 2m,„,e - ^— + V^'\r) , i/' ^ ^ __ . (2) 

For the IS state H', in the leading order, reads |l5lll6l| 

H' = Vna + Vnf + Vsi , (3) 

with the nonabelian and hyperfine terms given by 

CACpasifJ^y .. 87rC^a,(/i) 3 

^^^ - 2^^^^^^^ ' ^'f - 3ml ^ ^'> ' ^^^ 

and the spin-independent operator by 



Vsj = -^ + ^^^^ 5Hr) - ^^p- . (5) 



We shall now focus on the zeroth order Hamiltonian Ho. As mentioned the pole mass 
nipoie and singlet potential V^*^"-' in Hq are ill-defined inpQCD because of the renormalon 



divergences in the perturbative matching relations [17|, |1 
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FIG. 1: The improved potential V^^'{r) obtained by matching the BR potential and Cornell 
potential. 



where ^^(/i) denotes the strong coupling constant. 

The Hamiltonian in BR scheme can be obtained by substituting mpoie,V^^\r) with the 

Borel resummed quantities tubr, V^Qjl{r) which are defined by [A Isl. Il2j| 
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The Borel transforms in these equations have perturbative expansions about the origin. 
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where (3q denotes the one loop coefficient of the QCD j3 function, and about the renormalon 
singularity: 



\/(°)(6,l/r) 



Gri 



(1 - 26)1+'' 
(1 - 2hY+^ 



1 + ci(l - 26) + C2(l - 26)^ + . . . + (analytic part) , 

1 + ci(l - 26) + C2(l - 2hf + . . .1 + (analytic part) , (9) 



where v and Ci,C2 are exactly known constants. The renormalon residues G^, Cj, can be 
computed perturbatively following Refs. [^ ,£] . Combining the expansions Eq. (jH)) and Eq. 
(jH)) in bilocal expansions, a good description of the Borel transforms can be obtained from 
the known NNLO perturbations in Eq. ©. For details we refer the readers to |2l. I3lll2l|. 

Note that in the Borel summations (|7j) we have chosen separate scales, /i = m-g^ for 
the pole mass and /i = 1/r for the potential, which are, respectively, the optimal scales for 
the perturbative expansions (jH)). Since the renormalon ambiguities must be RG invariant 
we can choose the renormalization scales independently, without spoiling the cancellation 
of the renormalon ambiguities between the pole mass and the static potential. This is the 
crucial feature of our Borel resummation approach. 

Although the Borel resummed potential is expected to give a good description at short 
distance it will fail at long distance. To obtain the potential at all distances we follow the old 



idea of interpolating the pQCD potential to the confining linear potential. There are many 
implementations of this idea, from the original Richardson potential to the recent one based 
on renormalon cancellation [16|. Our implementation is as follows. We combine the BR 
potential valid at short distances with the Cornell potential at long distance by matching 
them at an intermediate distance tq by 

^BrVo) = Vbornell(ro) , ^V^BrVo) = ^Vbomell(ro) , (lO) 

where 

Vcornell(^) = « H V CT . (11) 

r 
Our choice of the Cornell potential as the long distance potential is based not only on its 
phenomenological success but its capability to fit the lattice data well [l9|. We take the 
linear slope as c = 0.18 GeV^, which is motivated by the lattice calculation [2u\. Note 
that the matching equations completely determine the two free parameters in the Cornell 
potential. Thus our hybrid singlet potential reads 

^(0) ^ f V^BR(r), r<ro, .^2) 

The potential obtained this way is continuous and smooth (see, Fig. 1) at the matching 
point Tq. 

The Borel resummed potential in pure QCD {rif = 0) agrees well with lattice calculations 
to distances as large as r ^ (700MeV)~^ [l^. Assuming that this behavior does not change 
drastically as the light quarks are introduced, we shall take the matching scale to be Tq^ = 1 
GeV. It will be shown that the dependence of the T(1S) binding energy on the matching 
scale is very small. 

As we shall see in the next section the incorporation of the long distance confining poten- 
tial is important even for such a spatially small system as T(1S), whose wave function turns 
out to have a sizable support where the difference between the confining potential and the 
perturbative potential is substantial. This implies that T(1S) is not completely perturbative 
even without the ultrasoft effects that will be discussed in Sec. IIVI 



III. SYSTEMATIC ERROR IN PURE PERTURBATIVE APPROACH 

Recently there have been higher order calculations of the T(1S) spectrum within pure 
pQCD framework (see for example, jl5|, |2l|, |22, |2^). In this approach the leading order 



binding energy can be obtained by solving the Schroedinger equation of the Hamiltonian, 
obtained by the perturbative expansion of 2mpoie + V^^\r) in Hq using the power expansions 
(jH)) at a fixed renormalization scale /i. Because of the renormalon cancellation between 2mpoie 
and V^'^''-' the perturbative expansion of the spectrum is free from renormalon divergence (note 
that, throughout the paper, divergences due to the sub leading renormalons are ignored). 
However, this approach can have a significant systematic error that cannot be seen from 
the convergence of the perturbative expansion or the renormalization scale dependence of 
the computed spectrum, since 2mpoie + V^^^ (r) , although converges perturbatively, cannot 
produce the confining potential however high order it has been calculated. 

To test whether the incorporation of the confining potential is important for T(1S) state, 
we calculate the ground state energy of Hq in BR as well as the perturbative schemes 




FIG. 2: The total energy in BR method (sohd hne), pQCD method (dashed hne) and zeroth order 
wave function r^^('^)(r)^ shifted by 8 for displaying purpose. 

TABLE I: The zeroth order total energy E^^'' by pure perturbative expansion method with differing set of 
scale /i, the unit is in GeV. 
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and compare them. We first plot 2mpoie + V^^''{r) in Fig. 2. in both BR scheme and 
in pQCD using NNLO calculations of the expansions ©• We put ^^g = 4.2 GeV and 
as{Mz) = 0.118 for both BR and pQCD schemes. In pQCD case the pole mass in the 
kinetic term were substituted by the NNLO pole mass in Eq. (jUj). The plot shows that the 
difference between these two schemes becomes substantial at r > lGeV~^, where the wave 
function for T(IS') has significant support, and because of this the computed spectrum 
is sensitive on the employed method. We list the spectrum in Table I computed by the 
pQCD method. Since the BR result is E^^^ = 9.574 GeV, we find the energy difference, 
which depends on the renormalization scale /i, is at least 100 MeV. This implies that the 
achievable accuracy in an approach based on pure pQCD with a fixed renormalization scale 
is at best 100 MeV, and that an incorporation of the confining potential at long distance is 
essential for precision calculation of T(1S) spectrum. 



IV. RECALCULATION OF THE LEUTWYLER-VOLOSHIN EFFECTS 



The binding energy of the color singlet quarkonium is also affected by the emission and 
absorption of the ultrasoft gluons, which occur in pNRQCD through the singlet-octet-gluon 
vertex. This effect is nonperturbative and difficult to calculate. However, when the typical 
time scale of the quarkonium, T ~ l/mv'^, is much smaller than the vacuum fluctuation scale, 
1/Aqcd, the quarkonium would 'see' the ultrasoft gluon as a constant random background 
field. The ultrasoft effect then becomes the Stark effect which in the case of quarkonium 
was first investigated by Voloshin 2^ and Leutwyler |2q . 

Whether this Stark effect approximation is applicable to T(1S) is debatable, but nonethe- 



less we shall assume in the following it is the case. As argued in Ref. [26[, even if this 
assumption fails, it is however likely that the Stark effect approximation will give an upper 
bound on the the ultrasoft effect, since the local condensate approximation gives an upper 
bound on the the two-point correlation function of the gluon field strength tensors. 

According to Leutwyler and Voloshin the Stark effect can be large on T(1S) state and 
this is known to give a large uncertainty to the spectrum obtained in pQCD approach. The 
calculations by Leutwyler and Voloshin employed Coulombic potentials for the singlet as 
well as octet channels. Obviously, these potentials must not be good representations of the 
true potentials, and thus the validity of the calculations based on Coulombic potentials is 
questionable. 

Our scheme of matching the Borel resummed potential at short distance with the confining 
potential at long distance is well-suited for the re-examination of this problem. In this 
section we recalculate the Stark effect using the hybrid singlet and octet potentials that were 
obtained by interpolations of the Borel resummed potential at short distance potentials and 
the confining potentials at long distance. 

Surprisingly, our result, which will be given shortly, shows that the Stark effect in T(1S) 
state is very small. This is fortunate since it removes the large uncertainty coming from the 
ultrasoft nonperturbative effect. 

The nonperturbative correction to the ground state energy is given by 

6E^P = ^{a,G%) J2{^,y^—-r^\^,) , (13) 

where \l/o, Eq are, respectively, the wave function and energy of the ground state of ifo, and 
Hq denotes the Hamiltonian for the octet channel: 

Ho = 2mBK-^ + Vo{r), (14) 

where Voir) denotes the static octet potential. 

With a recent preliminary calculation of the NNLO term J22| the octet potential is known 
perturbatively to NNLO. At short distance it is repulsive, but at long distance the octet 
potential is expected to have a confining linear potential. A hybrid octet potential that has 
good short distance behavior and a confining potential at long distance can be constructed 
in a similar fashion for the singlet potential. We first obtain the Borel resummed short 
distance potential employing the NNLO bilocal expansion, and then match it at the radius 
ro = IGeV"^ to the fitted potential V{r) = d+VSM3 - 0.289 r + 0.077 r2, which is obtained 
by fitting the lattice data for the octet potential [28|. The matching condition is that the 
potential be continuous at r = ro, and this determines the constant d. Note that the 
asymptotic linear slope of the fitted potential -\/0.077 ~ 0.277 GeV^ is roughly consistent 
with the bag model calculation, according to which the slope of the octet potential is given 



by a/7/4 times that of the singlet potential |22|. The octet potential obtained this way is 
plotted in Fig. 3. Now the nonperturbative correction Eq. (|T3|) can be easily evaluated by 
numerically solving the equation 

{Ho-Eo)\xl>=r'\^o> . (15) 
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FIG. 3: The octet potential, which is obtained by the Borel resummation at r < 1 GeV ^ and by 
fitting the lattice data at r > 1 GeV~^. 

V. RESULTS 



To obtain the binding energy we first solve the Schrodinger equation for Hq numerically 
to obtain the zeroth order binding energy E^^^ and the wave function \l/*^°'', and treat the 
interactions Vjyp, Vhf and Vsi as perturbations. 

Before performing the calculation, we have to deal with two problems caused by our 
choice of the inverse of the interquark distance as the renormalization scale for the potential. 
The hyperfine and spin-independent operators in Eqs. (jU Ej) contain the delta function 
S^{r), which will result in vanishing couplings as[l/{r —>■ 0)] in the numerator. In this 
case the inverse of the interquark distance as renormalization scale is clearly inappropriate. 
We solve this problem by noting that the renormalization scale for V^^\r) , V^^\r) can 
be independently chosen from the rest terms in H' since these two terms should be RG 
invariant, independently from the rest terms of H'. We therefore choose the renormalization 
scale /i = 1/r for V^^^ and Vna and a fixed value of /i for the hyperfine and spin-independent 
terms. The second problem is caused by the value of ^^(l/r) at the opposite condition, when 
r is large. The coupling ^^(l/r) obtained by the perturbative (3 function is not reliable when 
r is large. To avoid this difficulty, one can use a simple regulation |3C 



1^-^QCD — > a X tanh 
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qcd \ 






(16) 



in the running coupling as(l/r) = f^rAqcD), which is obtained using the four-loop /? 
function in MS scheme. The free parameter a should not be chosen too small a value so that 
the regularized coupling does not deviate too much from the perturbative coupling at short 
distances, nor too large a value to avoid an unphysical dip in the vicinity of the Landau 
pole. In our calculation, we set a = 0.55. We shall see that the computed T(1S) spectrum 
has only a mild dependence on this value. With this choice the regularized coupling changes 
little from the perturbative coupling when r is smaller than r = 1 GeV^^, and at large r the 
regulation suppresses the growth of the perturbative coupling and approaches to an infrared 
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FIG. 4: Regularized as{l/r). In infrared region as(l/r) tends to a fixed value. 

fixed value ^^(O) ^1.5 (see Fig. E}. This infrared behavior of the regularized coupling is 
interesting, considering that the effective running couplings in various dispersive approaches 
display freezing behavior in the infrared limit [^. Incidently, our q;s(0) ?« 1.5 is very close 
to the predicted infrared fixed value l//3o = 1.507 (at n/ = 4) of the effective coupling in 
"analytic perturbation theory" |3 2j| . 

Setting the input values of the strong coupling constant as{Mz) = 0.118, we solve the non- 
relativistic Schrodinger equation numerically with the hybrid potential, obtain the zeroth 
order binding energy E^'^^ and wave function \l/(°\ We treat the interactions Vna, Vhj and Vsi 
as perturbations and compute the perturbative corrections by calculating the expectation 
value 

5E = (vl/W I {Vna + Vhf + Vsi) \ v[/(°)). (17) 

The contribution from each term is given in Table |n] The values for the hyperfine and 
spin-independent contributions are with the renormalization scale /i = 3 GeV. 

The nonperturbative Leutwyler-Voloshin effect defined in Eq. ()13|) is computed by 
solving the equation (fT3j). Taking the gluon condensate (asGf^^,) = 0.02 ± 0.02 we find 
6E]\fP = 5.8 ± 5.8 MeV, which is almost an order of magnitude smaller than estimates based 

on Coulombic potentials, which range 30 100 MeV. We note that the smallness of this 

nonperturbative effect is very robust, depending little on the shape of the octet potential 
above the matching distance r^. We also report, without details, that Leutwyler-Voloshin 
effect on the hyperfine splitting is very small, at most 1.6 MeV. This result indicates that 
the ultrasoft nonperturbative effect, which is regarded a major difficulty in precision com- 
putation of quarkonium spectrum, is small in T(1S) state. 

Adding all these contributions as well as the charm mass effect, which shifts the binding 
energy about —20 ± 15 MeV [2], to ffi the measured mass of T(1S), we obtain the bottom 
quark BR mass, rriBR = 4.901 GeV, and the corresponding MS mass rriy^ = 4.211 GeV. 

It is noteworthy that the hyperfine and spin-independent contributions cancel largely. 
The sum of these two contributions is —20.5 + 30.4as(/^), in MeV, which shows that it has 
only a small renormalization scale dependence. Thus, our prediction of the 6-quark mass 
depends little on the renormalization scale chosen for the hyperfine and spin-independent 
terms. 

Our prediction of the hyperfine splitting can be found in Table IIIII and Table IIVI This 
result includes the leading as well as the subleading perturbative corrections for the hyperfine 
splitting [3^]. At the leading order the hyperfine splitting is quite sensitive on the value of 
the renormalization scale, but when the subleading correction is included it has a very 
small renormalization scale dependence for fi > 2 GeV. Taking /i = 3 GeV and as{Mz) = 
0.118, and the subleading correction as the theoretical uncertainty, we estimate the hyperfine 



10 



4.24 



4.22- 




0.114 0.116 0.118 \ 0.122 



FIG. 5: The extracted 6-quark MS mass vs as{Mz)- 

splitting to be 50 ± 8 MeV. 

In the following we discuss on the uncertainties on the extracted 6-quark mass in our 
scheme, (i) Uncertainty of as{Mz): With different input values of as{Mz), the extracted 
BR mass as well as MS mass will change. Fig. El shows the estimates of 6-quark MS mass 
against as{Mz). We can see from the figure, if the value of as{Mz) vary from 0.117 to 0.120, 
the MS mass will vary from 4.220 GeV to 4.192 GeV. (ii) Uncertainty of tq: Varying the 
matching point tq, at which the phenomenological potential is stitched to the perturbative 
BR potential, will cause some uncertainty. With about 10% variation of tq, the extracted b 
quark MS mass varies only within 3 MeV. This uncertainty can be completely ignored, (iii) 
Uncertainty of A: In our calculation, we chose the slope of the linear potential A = 0.18 GeV^. 
We shall assign 10% variation on A, which causes the MS 6-quark mass vary 4.211 ± 0.001. 
Again the error is small enough to be ignored. The T(1S) state is spatially small enough that 
its spectrum is not so sensitive on the shape of the confining potential, (iv) Uncertainty of 
the as regularization constant a: The regularization equation (fTB|) of a.^ involves the constant 
a, on which the infrared fixed value of ^^(l/r) depends strongly. For instance when we vary 
the value of a from 0.35 to 0.75, about 36% variation from the central value a = 0.55, the 
infrared fixed value of as{l/r) vary from 0.58 to 28.6, respectively. However, the variation 
of the extracted MS 6-quark mass is negligibly small, only from 4.205 to 4.217. This is 
similar to the small dependence on the slope A. Since only the long distance behavior of the 
regularized coupling is sensitive on a the spectrum of the T(1S) has only a small dependence 
on the value of a. 





TABLE II: Contributions to 15 state in unit of MeV. 




E(o) 


Vsi Vna Vhf 




9575 


-66 -113 58 


6 



VI. SUMMARY 

We investigated the T(1S) system using a hybrid Hamiltonian for heavy quarkonium, 
which was built by matching the Cornell potential at long distances to the Borel resummed 
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pQCD potential at short distances. This matching fixes completely the free parameters 
of the Cornell potential except for the slope of the linear term. The hybrid Hamiltonian 
employs an infrared sensitive, all-order resummed pole mass and static potential, and is 
free from the inevitable scale mixing problem present in approaches based on short- distance 
quantities. 

The incorporation of the long distance confining potential was shown to be essential for 
precision calculation of T(1S) system. The T(1S) is not a completely perturbative system, 
even without the ultrasoft contributions. We note that pure pQCD calculations of T(1S) 
spectrum at a fixed renormalization scale should contain a sizable systematic error of about 
100 MeV. 

The nonperturbative ultrasoft effects were reconsidered in the modeling of Leutwyler and 
Voloshin, employing the hybrid singlet and octet potentials. We find these effects are very 
small, contributing at most 6 MeV to the T(1S) spectrum, and negligible contribution to 
the hyperfine splitting which was found to be 50 ± 8 MeV at as{Mz) = 0.118. 

TABLE III: The hyperfine spHtting (leading order and next to leading order) vs as{Mz) at /i = 3 GeV, 
the unit is in MeV. 
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TABLE IV: The numerical results with different set of /i at as{Mz) ~ 0.118, the unit is in McV. 



^ 


1 


2 


3 


4 


5 


Vkf 


101 


69 


58 


52 


49 


^Vhf 


-81 


-21 


-8 


-3 


0.4 


Vhf + SVhf 


20 


48 


50 


50 


49 


Vsi 


-101 


-75 


-66 


-61 


-58 


"^MS 


4.207 


4.210 


4.211 


4.212 


4.212 
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